AD-A148  622  RESEARCH  IN  STOCHASTIC  PROCESSES(U)  NORTH  CAROLINA  UNIV  i/1 
AT  CHAPEL  HILL  DEPT  OF  STATISTICS  S  CAMBANIS  ET  HL. 

OCT  84  AFOSR-TR-84-1080  F49620-82-C-0009 


UNCLASSIFIED 


F/G  12/1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  •m*€AO  Of  STANOAAOS  -  t*A5  **  A 


UNCLASSIFIED  _ 


SECURITY  CLASSIFICATION  of  THIS  page 


report  security  classification 

UNCLASSIFIED 


2».  security  classification  authority 


ffl 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b  OECLASSIFICATION/OOWNGRADING  schedule 


a.  performing  organization  report  numberiSi 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERISI 

AFOSR-TR.  8  4- 


IBa.  NAME  OF  PERFORMING  ORGANIZATION  |6b.  OFFICE  SYMBOL  7*.  NAME  OF  MONITORING  ORGANIZATION 

University  of  North  w applicable i 

Carolina  Air  Force  Office  of  Scientific  Research 


6 c.  ADDRESS  iCity .  Slate  and  ZIP  Code t 

Statistics  Department 
Chapel  Hill  NC  27514 


Ba.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

AFOSR 


Sc  AOORESS  iCity.  Stale  and  ZIP  Code I 


7b.  ADDRESS  (City,  State  and  ZIP  Coda 

Directorate  of  Mathematical  &  Information 
Sciences,  Bolling  AFB  DC  20332-6448 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable  i 

m  F49620-82-C-0009 


10.  SOURCE  Of  FUNDING  NOS 


PROJECT  TASK  WORK  UNIT 

NO.  NO.  NO 


PROGRAM 

PROJECT 

TASK 

E  lEMENT  NO. 

NO. 

NO. 

Bolling  AFB  DC  20332-6448 

61102F 

2304 

A5 

11.  TITLE  (Include  Security  Classification/ 

RESEARCH  IN  STOCHASTIC  PROCESSES 


12.  PERSONAL  AUTHORlS)  _ _ _ ; _ 

S.  Cambanis.  R.J.  Carroll,  G.  Kallianpur,  and  M.R.  Leadbetter 


13*.  TYPE  OF  REPORT  I  13b.  TIME  COVERED  14.  DATE  OF  REPORT  Or .  .Mo  .  Z>«>7  15  PAGE  COUNT 

Interim  I  from  1/11/83  to  31/8/841  OCT  84  60 


16.  Supplementary  notation 


•>  'V  ■ 


17 

COSATl  CODES  | 

FIELD 

GROUP 

SUB  GR 

""  ” 

.. 

i.  c;s,  <  ■  ■  •  ■  t 


18-  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  6y  block  number / 


19.  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number/ 

■"^Research  was 'conducted^and  directed  in  the  area  of  stochastic  processes  by  three  of  the 
principal  investigators,  S.  Cambanis,  G.  Kallianpur,  and  M.R.  Leadbetter,  and  their 
associates,  and  in  statistical  estimation  and  inference, by  R.J.  Carroll  and  co-workers. 

A  summary  of  the  main  lines  of  activity  in  each  area  fc/flows  for  each  of  the  four 
principal  investigators.  More  detailed  descriptions  of  the  work  of  all  participants  is 
given  in  the  report.  |pae».^ 

Li  I  IU 

— CTE 


0112  FILE  COW 


DEC  1  4  1984 


20  DISTRIBUTION/AVAILABILITY  OF  abstract 
JNCLASSIFIEO/UNLIMITEO  £  SAME  AS  RPT  Z  OTIC  USERS  Z 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

MAJ  Brian  W.  Woodruff 


DO  FORM  1473,  83  APR 


22b  TELEPHONE  NUMBER 
i Include  \  rvo  C  ou* 

(202)  767-  L>0  ’7 


%  /.  ,%  , 


•*. 


0*F  CE  SVMBC 

n:-* 


CLASSIFIED 


EDITION  OF  1  JAN  73  IS  OBSOLETE 


HFOSR.TR.  84.10  80 


SCIENTIFIC  REPORT 

Air  Force  Office  of  Scientific  Research  Grant  AFOSR  F49620  82  C  0009 
Period:  I  November  1983  through  31  August  1984 


Title:  Research  In  Stochastic  Processes 


Co-Principal  Investigators: 


Professor  Stamatls  Cambanis 
Professor  Raymond  J.  Carroll 
Professor  Goplnath  Kalllanpur 
Professor  H.  Ross  Leadbetter 


Accession  For^ 

NTIS  GF.A&I 
DTTC  TAB 

Unannounced 
J-.r.t  if  lection. 


DiftrH  -.-.Mon/ 


B 


A  vc 

Dist 


Liability  Codes 
(Avail  and/or 
Special 


Approved  for  public  release  • 
distribution  unlimited. 


Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina  27514 


t  9  A  B* 


AIR  FORCE  OFF'.?:  OF  fCTENTIFIC  RESS.V^M  (  \FSC 
NOTICE  C?  "7'  TO  D?IC 

Thi3  t<  r  -  •  ,.,5 

tfpp  -  ■  .  ’  •  J-IZ. 

LlSt:  If  .  _ 

MATTHEi.  J. 

CONTENTS151 8f»  Technical  Inforaation  Division 


SUMMARY  OF  RESEARCH  ACTIVITY . 1 

RESEARCH  IN  STOCHASTIC  PROCESSES . 4 

Principal  Investigators: 

S.  Cambanls . 5 

G.  Kalllanpur . 10 

M.R.  Leadbetter . 15 

Visitors:  * 

O.A.  Bucklew . 18 

O.J.  Daley . 19 

R.A.  Davis . 20 

L.  de  Haan . 22 

T.  Funakl . 23 

W.N.  Hudson . 24 

R.  Jajte . 26 

R.L.  Karandlkar . 27 

A.F.  Karr . 29 

A.H.  Korezlloglu . 30 

H. H.  Kuo . 31 

H.O.  Nleml . 32 

M.  Pourahmadl . 33 

H.  Rootzen . 35 

J.  Roslnskl . 36 

D.  Tjpsthelm . 38 

RESEARCH  IN  STATISTICAL  ESTIMATION  AND  INFERENCE . 41 

Principal  Investigator: 

R.J.  Carroll . 42 

Visitor: 

H.  Schneider . 47 

PUBLICATIONS . 49 

CENTER  FOR  STOCHASTIC  PROCESSES  TECHNICAL  REPORTS . 52 

INSTITUTE  OF  STATISTICS  MIMEO  SERIES  TECHNICAL  REPORTS. .. .55 

STOCHASTIC  PROCESSES  SEMINARS . 56 

PROFESSIONAL  PERSONNEL . 58 

INTERACTIONS . 59 


SUMMARY  OF  RESEARCH  ACTIVITY 

Research  was  conducted  and  directed  In  the  area  of  stochastic  processes  by 
three  of  the  Principal  Investigators,  S.  Cambanis,  6.  Kalllanpur,  and  M.R.  Lead- 
better,  and  their  associates,  and  In  statistical  estimation  and  Inference  by 

R. J.  Carroll  and  co-workers.  A  summary  of  the  main  lines  of  activity  In  each 
area  follows  for  each  of  the  four  Principal  Investigators.  More  detailed  descrip¬ 
tions  of  the  work  of  all  participants  Is  given  In  the  main  body  of  the  report. 

STOCHASTIC  PROCESSES  * 

The  research  effort  In  stochastic  processes  was  a  major  part  of  a  substantial 
research  activity  organl2ed  as  the  Center  for  Stochastic  Processes  within  the 
Statistics  Department,  Involving  permanent  faculty,  visitors  and  students.  This 
organization  has  provided  the  framework  for  significant  Interaction  between  the 
participants— permanent  and  visiting..  In  addition  the  research  program  has  been 
enhanced  by  a  regular  seminar  series  (listed  by  speakers  later  In  the  report) 
which  has  provided  an  excellent  vehicle  for  exchange  of  current  research  Ideas. 

The  primary  means  for  dissemination  of  results  Is  the  Center's  Technical  Reports 
series,  containing  current  research  work  prior  to  formal  publication  In  journals.  To 
date,  80  technical  reports  leading  (or  expected  to  lead)  to  published  papers  have 
been  produced  by  the  participants.  Involving  research  results  In  a  wide  area  of 
stochastic  process  theory  and  applications.  The  main  areas  of  research  activity 
for  each  Principal  Investigator  and  co-workers  are  as  follows. 

S.  Cambanis.  Non-Gausslan  signals:  similarities  and  contrasts  between  Gaussian 
and  other  stable  signals;  Interpolation  of  harmonizable  stable  processes;  erqodic 
properties  of  stationary  stable  processes;  Wold  decompositions;  moment  inequalities 
for  stable  Integrals;  stochastic  Integral  representations  of  stable  processes 

with  paths  In  Banach  spaces;  stochastic  integrals  with  respect  to  Independent 


Increments  processes;  operator  stable  distributions.  Prediction  and  modeling 
of  time  series:  prediction  error  for  small  lags;  rate  of  convergence  of  finite 
linear  predictors;  autoregressive  representations  of  multivariate  series;  doubly 
stochastic  time  series  models;  estimation  In  nonlinear  time  series  models. 

Digital  processing  of  analog  signals:  sampling  designs  for  time  series;  designs 
for  estimating  random  Integrals  from  noisy  observations;  performance  of  discrete¬ 
time  predictors  of  contl nuous-tl me  processes. 

6.  Kalllanpur  Nonlinear  filtering  theory.  Stochastic  differential  equations. 
Finitely  additive  white  noise  theory:  Markov  and  robustness  properties;  non¬ 
linear  prediction  and  Interpolation.  Stochastic  differential  equations  In 
nuclear  spaces  and  application  to  neurophysiological  problems.  Stochastic  and 
Multiple  Wiener  Integrals.  Feynman  Integrals  and  abstract  Wiener  spaces.  Non- 
commutatlve  probability  theory.  Statistical  mechanics:  diffusions  approxima¬ 
tions  to  the  Boltzmann  equation.  Differential  geometry  in  statistical  problems 
of  stochastic  processes.  Brownian  functionals.  Stationary  processes  and 
prediction  theory. 

M.R.  Leadbetter  Extremal  theory:  extremes  and  local  dependence  in  stationary 
sequences;  extremes  In  dependent  sequences,  continuous  stochastic  processes, 
autoregressive  processes;  estimation  of  tails  of  distributions  from  extreme 
order  statistics.  Point  processes  of  exceedances:  compound  Poisson  limit  theorems 
for  high  level  exceedances,  point  processes  describing  the  joint  behavior  of 
order  statistics.  Estimation  for  point  processes:  maximum  likelihood  estimation 
In  the  multiplicative  Intensity  model.  Other  point  process  theory:  probability 
generating  functionals,  waiting  time  In  queues.  Estimation  and  structure  for 
stationary  processes:  probability  density  estimation. 


R.J.  Carroll:  Heteroscedastlclty:  Levene's  test  for  equality  of  variances, 
bounded-influence  estimation  In  heteroscedastlc  models,  diagnostics.  Binary 
regression:  errors-ln-varlables  models,  small  measurement  errors  asymptotics, 
robust  and  optimal  bounded- Influence  estimation  with  extensions  to  generalized 
linear  models.  Transformations  In  regression:  bounded  Influence  estimation, 
modified  likelihood  ratio  tests.  Accelerated  life  testing:  diagnostics  for 
Influential  points.  EM  algorithm  for  censored  data:  estimating  variance  In 
Buckley-James  model.  Sampling  plans:  numerical  analysis  of  plans  based  on 
prior  distributions  and  costs,  variable  sampling  plans  when  the  normal  distri¬ 
bution  Is  truncated.  Inventory  policies:  multi-item  Inventory  problem  using 
optimal  policy  surfaces. 


5 

STAMATIS  CAMBANIS 

The  work  briefly  described  here  was  developed  In  connection  with  problems 
arising  from  and  related  to  the  statistical  communication  theory  and  the 
analysis  of  stochastic  signals  and  systems,  and  falls  Into  the  following  three 
broad  categories: 

(I)  Non-Gaussian  signal  processing, 

(II)  Digital  processing  of  analog  signals. 

Items  1  to  3  below  belong  to  category  I.  Items  1  and  2  represent  continuing 
work  with  Drs.  Hardin  and  Heron.  Further  work  In  progress  on  non-Gausslan 
signals  will  be  reported  In  the  next  reporting  period;  some  jointly  with  Mr. 
Marques,  a  Ph.D.  student,  on  discrimination  between  stable  signals  and  detection 
of  sure  signals  In  stable  noise,  and  some  jointly  with  Dr.  Hudson  on  certain 
classes  of  Infinitely  divisible  signals.  Items  4  to  6  belong  to  category  II. 
Item  5  Is  joint  with  Dr.  Bucklew.  Item  6  represents  continuing  work  with  Dr. 
Masry  and  provided  the  Impetus  for  further  work  on  sampling  designs  which  will 
be  completed  and  reported  during  the  next  reporting  period. 

1.  Erqodlc  properties  of  stationary  stable  processes  [1] 

We  derive  spectral  necessary  and  sufficient  conditions  for  stationary 
symmetric  stable  processes  to  be  metrically  transitive,  mixing,  and  to  satisfy 
a  law  of  large  numbers,  and  we  show  how  to  estimate  their  covariation  functions. 
He  then  consider  some  Important  special  classes  of  stationary  stable  processes: 
Sub-Gausslan  stationary  processes  are  never  metrically  transitive,  and  we 
Identify  their  ergodlc  decomposition.  Stationary  stable  processes  with  a 
harmonic  spectral  representation  are  never  metrically  transitive.  In  sharp 
contrast  with  the  Gaussian  case,  and  some  Intriguing  facts  concerning  their 
ergodlc  decomposition  are  derived;  also  stable  processes  with  a  harmonic  spectral 
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representation  satisfy  a  strong  law  of  large  numbers  even  though  they  are  not 
generally  stationary.  For  doubly  stationary  stable  processes,  sufficient 
conditions  are  derived  for  metric  transitivity  and  mixing,  and  necessary  and 
sufficient  conditions  for  a  law  of  large  numbers. 

2.  Wold  decompositions  of  symmetric  stable  sequences  [2]. 

This  paper  Introduces  a  notion  of  Wold,  or  orthogonal,  decomposition  of 
non-second-order  processes  and  then  gives  necessary  and  sufficient  conditions 
for  Its  existence  In  the  case  of  symmetric  stable  sequences.  The  orthogonality 
used  here  (the  usual  orthogonality  being  Inapplicable  In  the  non-second-order 
case)  Is  known  as  James'  orthogonality  —  a  non-synvnetric  notion  defined  for 
Banach  spaces  which  Is  In  many  ways  the  most  appropriate  for  our  situation. 

We  show  there  Is  an  Intimate  relation  between  orthogonality  of  symmetric 
stable  variables.  This  Is  exploited  to  show  that  a  stable  sequence  has  a  Wold 
decomposition  If  and  only  If  the  regression  of  each  random  variable  on  the 
preceding  ones  Is  linear.  Also  given  are  conditions  for  the  coincidence  of  this 
decomposition  with  the  "non-linear"  Wold  decomposition  (which  Is  shown  to  always 
exist),  and  an  equivalent  spectral  condition  for  Independence  In  the  Wold 
decomposition. 


3.  Similarities  and  contrasts  between  Gaussian  and  other  stable  signals  [3], 
Gaussian  stochastic  signals  are  the  simplest  and  most  thoroughly  studied 
among  all  stable  stochastic  processes,  but  are  they  typical  or  atypical?  Viewed 
another  way,  are  signal  processing  techniques  designed  for  Gaussian  signals 
robust  with  respect  to  mild  departures  from  normality  within  the  Important  class 
of  stable  stochastic  signals?  Different  answers  to  this  question  are  given  in 
a  variety  of  contexts  such  as  regression,  linear  estimation,  prediction  and 


filtering,  parameter  estimation,  signal  detection,  and  system  Identification. 
Sharp  contrasts  as  well  as  certain  similarities  between  Gaussian  and  other 
stable  processes  emerge.  Before  considering  specific  problems,  stable  stochastic 
signals  are  Introduced  and  discussed,  with  particular  attention  to  stationary 
stable  signals.  The  aim  Is  to  survey  recent  results  and  to  focus  attention  on 
open  problems. 

4.  Sampling  designs  for  time  series  [4]. 

In  practice  a  time  series  Is  observed  only  at  a  finite  number  of  points 
(the  sampling  design)  and  based  on  these  observations  an  estimate  or  a  statistic 
is  formed  for  use  In  the  problem  at  hand.  The  statistician  may  be  free  to  choose 
the  sampling  points,  or  part  of  the  sampling  mechanism  may  be  Imposed  on  the 
statistician  who  then  controls  only  certain  parameters;  e.g.  periodic  sampling 
Is  Imposed  where  the  period  Is  controlled  by  the  statistician,  or  sampling  at  the 
times  of  occurrence  of  a  Poisson  stream  of  events  is  Imposed  but  the  statistician 
has  control  over  Its  rate.  How  can  the  statistician  choose  the  best  design  of 
a  given  sample  size,  or  determine  the  sample  size  of  a  certain  kind  of  deslan 
required  to  achieve  a  given  performance?  These  questions  are  considered  in  the 
context  of  three  specific  problems  involving  time  series:  the  estimation  of  a 
weighted  average  of  a  random  quantity,  estimation  of  regression  coefficients,  and 
the  detection  of  signals  In  noise. 

The  set  up  here  differs  In  two  Important  ways  from  the  classical  set  up. 

All  observations  are  taken  from  a  fixed  region  and  so,  especially  for  large 
sample  sizes.  It  Is  not  realistic  to  assume  lack  of  correlation;  hence,  observa¬ 
tions  form  a  correlated  time  series.  Also  repeated  sampling  at  the  same  point 
Is  not  allowed,  and  only  one  realization  of  the  time  series  is  available. 

We  consider  both  deterministic  and  random  sampling  designs,  where  either 
optimal  estimators  and  sufficient  statistics  are  employed,  or  much  simpler 


estimators  and  statistics  are  employed  Instead.  Finding  optimal  designs  of  a 
given  sample  size  turns  out  to  be  a  very  difficult  task,  which  can  be  accomplished 
only  for  certain  specific  covariance  structures  or  certain  simple  sampling 
designs,  such  as  simple  random  sampling.  Finding  sampling  designs  which  for 
large  sample  size  perform  like  the  best  designs,  is  an  easier  task,  which  can 
be  accomplished  for  broad  classes  of  covariance  structures  and  different  designs, 
such  as  unconstrained,  median  and  stratified,  the  latter  two  in  fact  using  the 
simpler  kind  of  estimators  and  statistics.  A  unified  approach  is  presented  here 
and  the  results  on  the  simpler  estimators  of  regression  coefficients  are  new. 


5.  Sampling  designs  for  estimating  random  integrals  with  observation  noise.  [5]. 

We  consider  the  problem  of  estimating  a  weighted  average  of  a  random 
quantity  using  observations  corrupted  by  noise  at  a  finite  number  of  sampling 
points,  and  we  study  the  effect  of  the  observation  noise  on  the  sampling  design 
and  on  Its  performance.  The  cases  of  additive  Independent  noise,  of  nonlinear 
distortion  plus  noise,  and  of  quantization  noise  are  considered  in  detail. 

6.  Performance  of  discrete-time  predictors  of  continuous-time  processes  [6], 

He  study  the  asymptotic  performance  of  linear  predictors  of  continuous- time 
stationary  processes  from  observations  at  n  sampling  instants  on  a  fixed  observa¬ 
tion  Interval.  He  consider  both  optimal  and  simpler  choices  of  predictor 
coefficients;  and  uniform  sampling,  as  well  as  nonuniform  sampling  tailored 
to  the  statistics  of  the  process  under  prediction.  We  concentrate  on  stationary 
processes  with  rational  spectral  densities  and  numerical  examples,  depicting 
small  sample  size  performance  In  addition  to  asymptotics,  are  given  for  cases 
with  no  and  with  one  quadratic  mean  derivative. 
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GOPINATH  KALLIANPUR 

1.  Finitely  additive  approach  to  nonlinear  filtering  (with  R.L.  Karandikar)[T] 


This  Is  the  continuation  of  a  project,  begun  In  1982,  of  developing  a 
completely  new  approach  to  nonlinear  filtering  theory  when  the  signal  and  noise 
are  Independent.  Significant  progress  was  achieved  during  the  period  of  the 
present  report  on  the  following  problems: 

(I)  The  optimal  filter  In  the  general  white  noise  model  Is  shown  to  be  a 
Markov  process.  More  precisely.  It  Is  shown  that  Ft(y)  and  rt(y)  -  the 
normalized  and  unnormalized  conditional  distribution  (where  y  is  the  white 
noise  observation)  regards  as  measure-valued  processes  either  on  the  quasl- 
cyllndrlcal  probability  space  (E,E,a)  or  as  defined  on  (H,C,n)  are  Markov 
processes.  (For  definitions  and  details,  see  Technical  Report  No.  44.  [2]) 

The  principal  difficulty  until  now  has  been  the  lack  of  a  suitable  trans¬ 
formation  formula.  The  question  is  solved  by  paying  more  attention  to  the 
underlying  theory,  by  discovering  a  more  versatile  definition  of  conditional 
expectation  (See  Section  2). 

(II)  Robustness  results:  A  crucial  achievement  of  the  white  noise  theory  Is 
Its  robustness  properties.  Some  of  the  results  were  reported  In 

[8]  but  the  best  were  proved  with  the  aid  of  the  stronger  theory 

worked  out  in  [3]  and  are  Included  there.  In  addition  to  robustness, 

the  consistency  of  the  finitely  additive  filtering  approach  with  the  conventional 

theory  based  on  martingale  calculus  and  Ito  stochastic  differential  equations 

have  been  established.  A  typical  result  may  be  briefly  stated  as  follows: 

Let  y  *  h  (X.)  +  e  ,  0  <  s  <  T  be  the  finitely  additive  filtering  model 

defined  on  the  quasl-cyllndrlcal  probability  space  (E,E,ot)  *  (n,A,H)  ®  (H,C,m) 

2  d 

where  e  Is  Gaussian  white  noise,  H  *  L  [0,T]  and  (X$)  Is  the  1R  -valued 


.1 
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signal  process  defined  on  (n,A,n).  The  stochastic  calculus  model  Is  Y+  * 

it  1 

■  I  hs(Xs)ds  ♦  Wt,  («)  being  a  standard  Wiener  process.  Let  pt(x,n)  be 
the  unnormallzed  conditional  density  of  given  y,  obtained  as  the  unique 
solution  of  the  finitely  additive  Zakal  equation.  Then  (a)  there  exists  a 
continuous  function 

P  (•  »Z):  C0([0,T];»'n)  -  C([0,T]  *  lRd) 
such  that  for  all  n  c  H, 

A  d  ft 

Pt(x,Z)  «  pt(x,n)  .  0<t<T,xcR  where  Zt  =  j  ^ds*  ; 

a 

(b)  Pt(x,Y)  Is  a  version  of  the  unnormallzed  conditional  density  of  Xt 

V 

are  given  F^. 

The  precise  conditions  under  which  the  above  theorem  Is  proved  are  given  In  [2]. 
The  result  shows  that  the  robustness  results  of  the  stochastic  calculus  approach 
(e.g.  the  results  of  Pardoux)  can  be  obtained  more  simply  from  the  white  noise 
theory  under  slightly  weaker  conditions. 

2.  White  noise  calculus  (with  R.L.  Karandikar) 

A  full-fledged  development  of  the  theory  behind  the  new  approach  to 
nonlinear  filtering  theory  (described  In  Sec.  1)  was  undertaken  In  [3]  which  Is 
a  first  draft  of  a  forthcoming  monograph.  The  authors  Introduce  Integration 
with  respect  to  a  quasi -cylinder  probability  measure  (QCP)  of  suitable  classes 
of  cylinder  functions,  followed  by  a  definition  of  absolute  continuity  of 
QCP's.  Although  the  notions  of  a  representation  and  lifting  associated  with 
a  QCP  are  similar  to  Ideas  of  I.E.  Segal  and  L.  Gross,  the  problems  of 
nonlinear  filtering  are  sufficiently  different  as  to  require  a  detailed  and 
Independent  treatment.  Another  most  useful  and  Important  concept  Is  that  of  a 
quasi -cylindrical  mapping  (QCM)  which  has  most  of  the  desirable  properties  of 
a  random  variable  In  the  countably  additive  theory.  A  QCM  defined  on  a 
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quasi  cylindrical  probability  space  induces  a  QCP  on  the  range  space  and  an 
associated  Induced  representation.  Furthermore  it  yields  a  definition  of 
conditional  expectation  more  inclusive  and  useful  than  the  one  to  be  found 
in  an  earlier  paper  of  ours. 

3.  Feynman  integrals  (with  D.  Kannan  and  R.L.  Karandikar)  [4] 

Reported  under  Karandikar *s  work. 

4.  Application  of  differential  geometry  in  statistical  problems  of  stochastic 
processes  (with  H.H.  Kuo)  [5] 

Reported  under  Kuo's  work. 

5.  Regularity  property  of  Donsker's  delta  functional  (with  H.H.  Kuo)  [6] 

Reported  under  Kuo's  work. 

6.  Stochastic  differential  equation  models  for  neuronal  behavior  (with  R. 


Wolpert,  Duke  University).  [7] 

This  Is  a  continuation  of  earlier  work  (described  in  the  last  report)  in 
which  changes  in  the  voltage  potential  of  a  spatially  distributed  neuron  was 
treated  as  a  stochastic  process  taking  values  in  a  suitable  nuclear  space.  A 
nuclear  space-valued  Ornsteln-Uhlenbeck  process  was  obtained  as  a  diffusion 
approximation.  In  the  present  work,  see  [7]  more  realistic  models  are  considered 
(under  the  mathematically  simplifying  assumption  that  the  spatial  extension  of 
the  neuron  is  neglected)  which  take  Into  account  such  nonlinear  phenomena  as 
reversal  potentials,  or  the  situation  when  the  sizes  of  the  postsynaptic 
Impulses  depend  on  the  state  of  depolarization  of  the  membrane  potential. 
Non-gausslan  diffusion  approximations  to  neuronal  behavior  are  obtained.  As  a 
special  case  we  obtain  the  diffusion  approximation  to  the  Wan-Tuckwell  equation 
with  reversal  potentials.  The  latter  problem  leads  to  an  interesting  investigation 
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of  boundary  behavior. 

The  next  phase  of  the  work  Is  to  extend  the  above  theory  to  spatially  dls 
trlbuted  neurons.  Also  on  the  agenda  Is  a  satisfactory  description,  from 
the  stochastic  point  of  view,  of  the  celebrated  Hodgkin-Huxley  equations. 

Ph.D.  students  under  Goplnath  Kalllanpur 

Victor  Perez-Abreu:  Research  on  his  Ph.O.  thesis  Is  on  the  following  topics: 

(1)  Product  stochastic  measures:  A  new  approach  to  product  stochastic 

measures  has  been  developed  which  systematically  exploits  the  use  of  symmetric 
tensor  product  Hilbert  spaces.  Comparison  with  Engel's  work  Is  discussed  and 
other  examples  are  considered. 

(2)  Multiple  Wiener  Integrals  for  nuclear  space  valued-processes: 

Multiple  Wiener  Integral  expansions  and  stochastic  Integral  representa¬ 
tions  have  been  obtained. 

Soren  Christensen:  Stochastic  Differential  Equations  for  nuclear  space-valued 
processes  and  applications.  Research  on  his  Ph.D.  thesis  Is  on  extensions  and 
generalizations  of  Kalllanpur  and  Wolpert's  work  on  stochastic  differential 
equation  models  for  spatially  distributed  neurons.  The  case  of  a  semigroup 
(T^)  which  Is  not  necessarily  Hllbert-Schmldt  Is  considered.  Weak  convergence 
to  diffusion  approximations  (under  various  conditions)  Is  being  Investigated. 

Hans  Hucke:  Nonlinear  prediction  and  Interpolation  problems:  Hucke  Is 
extending  the  (finitely  additive)  white  noise  theory  of  Kalllanpur  and  Karandl 
kar  with  a  view  of  applying  it  to  nonlinear  prediction  and  smoothing  problems. 
General  signal  processes  will  be  considered.  He  will  also  Investigate  robust¬ 
ness  properties  for  these  problems  and  hopes  to  obtain  results  similar  to 
those  derived  for  nonlinear  filtering  theory. 


V  s*  V  O  «l»  O  O 
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M.R.  LEADBETTER 

During  this  reporting  period,  research  was  continued  or  Initiated  In  the 
following  main  areas:  (1)  continued  effort  In  the  theory  of  point  processes 
associated  with  extremal  theory  for  stochastic  sequences  (2)  extremal  theory 
for  non-statlonary  sequences,  (3)  further  development  of  extremal  theory  In 
continuous  time  and  (4)  function  estimation  for  stationary  processes.  The 
effort  In  each  is  described  below. 

(1)  Point  processes  associated  with  extremal  theory. 

Work  continued  with  J.  Husler  and  T.  Hslng  on  the  detailed  structure  of  the 
point  process  of  high  level  upcrosslngs,  extending  the  results  previously 
reported  In  [1],  Two  different  approaches  give  related  results.  One  of  these, 
described  In  the  report  [2],  gives  sufficient  conditions  for  convergence  of  the 
exceedance  point  processes  to  a  compound  Poisson  process.  The  other,  involving 
a  direct  "Laplace  Transform"  approach,  yields  a  complete  characterization  of  the 
limiting  point  process,  under  natural  conditions.  This  work  Is  to  be  reported 
In  [3]. 

Extensive  work  has  been  done  also  Into  so  called  "complete  convergence" 
results  Involving  the  point  process  consisting  of  the  values  of  a  stochastic 
sequence  plotted  (with  appropriate  normalizations)  In  the  plane.  Such  results 
give  a  comprehensive  picture  of  the  extremal  properties  of  the  sequence  -  including 
the  joint  asymptotic  distributional  behavior  of  the  maximum  and  other  extreme 
order  statistics.  A  theory  generalizing  that  of  Mori  [4]  has  been  obtained,  in 
which  the  possible  limiting  point  processes  are  fully  characterized.  This  has 
been  partially  reported  In  [5],  and  further  report  and  journal  publication  Is 
planned. 


(2)  Extremal  theory  for  non-statlonary  sequences 

Presently  available  extremal  results  for  non-statlonary  sequences  have 
typically  been  obtained  by  generalizing  those  for  stationary  cases  in  relatively 
restricted  and  obvious  ways.  An  Investigation  has  begun  taking  a  more  systematic 
approach,  first  considering  classes  of  normal  non-statlonary  processes  with 
rather  general  correlation  structure.  It  Is  planned  that  some  of  these 
results  will  be  Included  In  [6]  together  with  results  of  H.  Rootzen  on 
Markov  sequences. 

(3)  Extremal  theory  in  continuous  time. 

A  rather  satisfying  theory  for  extremes  of  continuous  parameter  stationary 
processes  has  been  developed  In  recent  years  and  reported  e.g.  in  the  book  [7]. 
Some  attention  has  been  given  In  the  current  reporting  period  towards  extending 
this  theory  to  cases  of  higher  local  dependence.  In  particular  this  involves 
the  attempt  to  parallel  the  discrete  time  theory  of  high  level  exceedances 
(cf.  (1)  above)  for  high  level  upcrosslngs  In  continuous  time.  However,  there 
are  also  features  which  appear  In  continuous  time  which  have  no  Interesting 
analogue  In  the  discrete  context  (due  to  the  presence  of  arbitrarily  high 
correlation  at  close  time  points).  These  are  being  investigated  -  considerable 
work  remains  to  be  done  to  obtain  definitive  results. 

(4)  Function  estimation  for  stationary  processes. 

The  research  on  this  topic  described  In  the  previous  report  has  been  further 
developed  In  the  current  period.  This  has  resulted  in  a  fairly  comprehensive 
account  of  consistency  and  asymptotic  distributional  properties  of  smoothed 
estimates  of  the  probability  density  of  a  stationary  sequence  or  continuous 
parameter  process.  This  work  (with  J.  Castellana)  has  been  reported  In  [8]  and 
submitted  to  0.  Stoch.  Proc.  for  publication. 
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JAMES  A.  BUCKLEW 


Dr.  Bucklew  completed  the  work  described  In  the  following  reports. 


1.  The  prediction  error  for  small  time  lags  Into  the  future  [1]. 

Explicit  expressions  are  derived  for  the  derivatives  at  zero  lag  of 
the  mean  square  prediction  error  of  a  class  of  stationary  processes,  which 
Includes  those  with  rational  spectral  density.  A  sufficient  condition  on 
the  spectral  density  Is  given  for  a  stationary  process  to  belong  to  this 
class,  and  It  Is  shown  that  the  stationary  process  with  triangular  covariance 
does  not  belong  to  the  class. 

2.  Sampling  designs  for  estimating  random  Integrals  with  observation  noise  [2]. 

We  consider  the  problem  of  estimating  a  weighted  average  of  a  random 
quantity  using  observations  corrupted  by  noise  at  a  finite  number  of  sampling 
points,  and  we  study  the  effect  of  the  observation  noise  on  the  sampling 
design  and  on  Its  performance.  The  cases  of  additive  independent  noise,  of 
nonlinear  distortion  plus  noise,  and  of  quantization  noise  are  considered  In 
detail . 
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Dr.  Daley  conducted  research  In  the  theory  of  queues  and  point  processes. 
This  work  Is  reported  in  [1]  and  [2],  whose  abstracts  are  as  follows: 

1 .  A  lower  bound  for  mean  characteristics  in  E^/G/1  and  GI/E^/l  queues  [1 ] . 

The  mean  stationary  waiting  time  EW  of  a  stable  E^/G/l  queuing  system 
depends  on  the  service  distribution  through  Its  first  two  moments  and  k-1  roots 
of  an  equation  Involving  Its  transform.  Vainshtein's  (1983)  work  on  delimiting 
the  region  where  such  roots  lie  Is  Improved  and  leads  to  a  bound  on  EU  that  Is 
always  tighter  than  other  known  lower  bounds.  The  methods  can  also  be  applied 
to  GI/E^/1  systems,  but  the  bounds  are  not  as  tight. 

2.  On  the  probability  generating  functional  for  point  processes  [2]. 

An  extended  probability  generating  functional  (p.g.fl.)  G  Is  defined 
for  a  point  process  N  on  a  complete  separable  metric  space.  The  distribution 
of  N  Is  uniquely  determined  by  G  over  a  small  class  of  functions.  Continuity 
results  are  reviewed  and  used  In  furnishing  a  complete  p.g.fl.  proof  of  the 
mixing  properties  of  certain  stationary  cluster  processes. 
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RICHARD  A.  DAVIS 

Dr.  Davis,  In  joint  work  with  Resnick,  Investigated  the  problem  of  estimating 
the  tall  of  a  distribution  function  using  extreme  order  statistics.  The 
aoDroach  assumes  one  knows  only  that  the  unknown  distribution  is  in  some 
domain  of  attraction  of  an  extreme  value  distribution.  The  results  of  this 
study  are  contained  In  the  paper  [1]  which  will  appear  December  1984  in  the 
Annals  of  statistics. 

Davis  and  Resnick  considered  the  large  sample  behavior  of  the  sample 

00 

covariance  and  correlation  functions  of  the  moving  average  X.  =  J  c.Zt  . 

1  j=-»  J  3 

where  (Zt>  Is  lid  with  coimion  distribution  In  the  domain  of  attraction  of  a 
stable  law  with  index  o,  0  <  a  <  2.  If  0  <  a  <  2,  E| Z^ |a  <  ®  and  the  distribu¬ 
tion  of  | Z-|  ]  and  IZ^I  are  tall  equivalent,  then  the  sample  correlation  function 
of  {Xt>  suitably  normalized  converges  In  distribution  to  the  ratio  of  two 
dependent  stable  random  variables.  On  the  other  hand  If  E|Z^ |tt  *  ®,  the  limit 
distribution  Is  the  ratio  of  two  Independent  stable  variables.  The  derivations 
of  these  results  rely  heavily  on  point  process  techniques.  In  the  a  =  2  case,  the 
sample  correlations  are  shown  to  be  asymptotically  normal  which  extends  the 
classical  result.  These  results  are  written  up  In  the  papers  [2]  and  [3]. 

Davis  and  Yao  considered  the  problem  of  testing  a  shift  in  mean  in  a  sequence 
of  Independent  normal  random  variables  In  the  paper  [4].  The  asymptotic  operating 
characteristics  of  the  likelihood  ratio  test  In  testing  HQ  (no  change  point) 
vs.  Hj  (change  point)  were  derived.  It  was  shown  that  the  asymptotic  null 
distribution  Is  closely  related  to  the  extreme  value  of  the  Ornstein-Uhlenbeck 
process  and  converges  In  distribution  to  the  double  exponential  extreme  value 
distribution. 
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LAURENS  de  HAAN 

Or.  de  Haan  conducted  research  In  extremal  and  related  theory,  with 
particular  reference  to  relationships  between  the  operations  of  "maximum" 
and  "sum".  This  work  Is  reported  In  [1]  which  discusses  stationary  processes 
that  are  max-autoregressive  In  one  time  direction  and  sum-autoregressive  In 
the  other. 
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'  Diffusion  approximations  to  the  Boltzmann  equation  [1.2]. 

Dr.  Funakl  has  continued  his  work  on  problems  of  statistical  mechanics, 
specifically,  the  study  of  the  Boltzmann  equation.  A  jump  process  Is  associated 
with  the  spatially  homogeneous  Boltzmann  equation  by  solving  a  martingale 
problem.  With  suitable  scaling,  this  process  is  shown  to  converge  weakly  to 
a  stochastic  process  related  to  a  nonlinear  diffusion  equation  known  as  the 
Landau  equation. 

In  a  similar  context,  existence  and  uniqueness  theorems  have  been  obtained 
for  the  martingale  problem  associating  diffusion  processes  with  a  type  of 
nonlinear  parabolic  equation. 
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WILLIAM  N.  HUDSON 


■ 

Dr.  Hudson  continued  his  work  on  multivariate  operator  stable  distribu¬ 
tions  and  Initiated  the  study  of  stochastic  Integrals  with  respect  to  general 
Independent  Increments  processes. 

1.  The  symmetry  group  and  exponents  of  operator  stable  probability 
distributions  [1]. 

It  Is  shown  that  there  exist  exponents  of  an  operator  stable  distribution 
which  commute  with  every  operator  In  the  distribution's  symmetry  group.  These 
exponents,  together  with  a  new  norm,  lead  to  useful  simplifications  in  the 
representation  of  the  distribution's  Levy  measure. 

3.  Stochastic  Integrals  with  Respect  to  Independent  Increment  Processes  [2]. 

The  stochastic  Integral  /^V(s)d?(s)  Is  defined,  where  £  Is  a  stochastically 
continuous  process  with  independent,  but  not  necessary  stationary.  Increments, 
without  a  Gaussian  component,  and  which  Is  adapted  to  a  filtration,  and  V  is  a 
stochastic  process  which  is  adapted  to  the  same  filtration.  If  the  paths  of 
C  are  not  of  bounded  variation  over  [0,t],  the  Integrands  V  are  assumed  to  be 
predictable  and  to  satisfy  the  condition 

/  {(xV(s))2  a  |xV(s)|}dM(x,s)  <  »  a.s. 

[-a,a>[0,t]- 

* 

for  some  a  >  0  where  M  denotes  the  jump-time  Levy  measure  of  5.  This  class 
is  larger  than  that  produced  by  the  usual  martingale  approach  and  was  first  used 
by  Kallenberg  [4]  who  assumed  s  had  stationary  Increments. 

The  results  concerning  approximation  of  Integrands  by  simple  processes 
and  the  existence  of  the  Integral  contain  new  features.  An  elementary  method 
is  used  to  obtain  simple  processes  which  approximate  the  Integrands.  An  argu¬ 
ment  of  Kallenberg  Is  extended  to  show  the  existence  of  this  Integral.  This 


stochastic  Integral  Is  also  shown  to  be  Indistinguishable  from  the  pathwise 

Lebesgue-Stleltjes  integral  whenever  both  Integrals  exist.  The  Integral  is 

shown  to  possess  the  usual  properties  and  a  complex- valued  exponential 

martlnnale  is  associated  with  the  integral. 

Finally,  attention  Is  focused  to  nonrandom  Integrands,  and  a  formula  is 

obtained  for  their  characteristic  function.  The  class  of  integrands  turns  out 

to  be  as  large  as  possible  consistent  with  this  formula,  which  provides  an 

easy  way  to  show  that  such  Integrals  are  themselves  stochastic  processes  with 

* 

Independent  increments.  The  jump-time  Levy  measure  Is  easily  obtained  and  used 

to  show  that  If  V1  and  V2  are  two  integrate  processes,  then  /J  V-,<fc  and 

/q  VgdCare  independent  random  variables  If  and  only  if  (sJVgtsJ^O 

v-a.e.  where  u(A)  *  /  (x2  a  l)dM(x,s). 

RxA 
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RYSZARO  JAJTE 


Strong  limit  theorems  In  von  Neumann  algebras.  [1] 

Several  classical  polntwlse  convergence  theorems  have  been  extended  to 
the  context  of  von  Neumann  algebras.  The  dynamics  of  a  quantum  mechanical 
system  S  (described  by  a  von  Neumann  algebra  A  and  a  state  <f>)  are  given  by  a 
one  parameter  group  (a^)  of  automorphisms  of  A.  It  Is  of  great  Interest  to 
know  the  asymptotic  behavior  of  the  ergodlc  averages 

(*)  ^  |  o  “kM  (ne1.2...)  asn^« 

for  x  «  A.  Theorems  related  to  this  problem  are  proved. 
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RAJEEVA  L.  KARANDIKAR 


1.  Finitely  additive  approach  to  nonlinear  filtering  theory  [1]  (with  G. 
Kalllanpur):  Reported  under  Kalllanpur's  work. 

2.  White  noise  calculus  [2]  (with  G.  Kalllanpur):  Reported  under  Kalllanpur's 
work. 

3.  Feynman  Integrals  [3]  (with  G.  Kalllanpur  and  D.  Kannan). 

Continuing  earlier  work  on  this  subject,  a  unified  approach  is  presented: 
Sequential  Feynman  integrals  are  defined  for  classes  of  functions  on  Hilbert 
space  and  on  an  abstract  Wiener  space.  A  Cameron-Martin  formula  is  proved  for 
analytic  and  sequential  Feynman  Integrals  for  the  classes  Gq(B)  and  Gq(H).  A 
function  g  on  H  belongs  to  Gq(H)  If  It  is  of  the  form 

g(h)  «  eV2(h,Ah)  j  (k»h)<ju(k)  where  q  (fO)  Is  a  real  number, 

li  Is  a  complex,  bounded  Borel  measure  on  H  and  A  Is  a  self-adjoint,  trace  class 
operator  on  H  such  that  I  +  ^  A  has  a  bounded  Inverse.  The  class  Gq(B)  has  an 
analogous  definition. 

4.  limiting  distributions  of  functionals  of  Markov  chains,  (with  V.G.  Kulkarni) 

If  Xn  Is  a  Markov  chain  with  stationary  transition  probabilities  and 

*  f(X  .....X^)  (k  then  depends  on  Xn  In  a  stationary  way  in  the 

sense  that  P(Y„  «  AlXj  does  not  depend  on  n.  Sufficient  conditions  for  Y„ 
n  1  n  n 

to  have  a  limiting  distribution  are  derived  for  the  following  two  cases: 

(1)  {Xn,  n  ^  0}  has  a  limiting  distribution,  (11)  {Xn,  n  >_ 0}  does  not  have  a 
limiting  distribution  and  exits  every  finite  set  with  probability  one.  Several 
examples  are  considered  Including  that  of  a  nonhomogeneous  Poisson  process  with 
periodic  rate  function  for  which  the  limiting  distribution  of  the  Inter-event 
times  Is  obtained. 
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Dr.  Karr  pursued  various  research  topics  In  the  general  area 
of  statistical  Inference  for  stochastic  point  processes. 

He  completed  the  technical  report  [1],  which  has  been  accepted  for  publica¬ 
tion  In  the  Anna! s  of  Statistics.  Further  work  dealt  with  strong  approximation 
of  Poisson  processes.  Inference  for  Poisson  processes  with  one  observable 
comment  and  Inference  for  point  processes  given  integral  data.  Karr  also 
completed  Chapters  3-8  of  his  forthcoming  book  Point  Processes  and  their 
Statistical  Inference. 
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AHMET  HAYRI  KOREZUOGLU 

1 .  Stochastic  Integration  for  operator-valued  processes  on  Hilbert  spaces  and 
on  nuclear  spaces.  [1]  (with  C.  Marti  as). 

The  aim  of  the  work  Is  to  give  a  representation  of  distribution-valued 
square-lntegrable  martingales  by  means  of  stochastic  Integrals  with  respect 
to  a  given  distribution-valued  martingale.  The  problem  of  representation  of 
Hllbert-space  valued  martingales  Is  reconsidered  so  as  to  Include  stochastic 
Integration  of  Hilbert- Schmidt  operator  valued  processes. 

2.  An  approximation  for  the  Zakai  equation  by  a  finite  difference  equation  [2]. 

Consider  the  filtering  problem, 

*t  -  4  £  j<vds  4  £  *>(xs>dBs  • 

yt  -  £  c(Xs)ds  ♦  Wt 

where  Xq,  B  and  U  are  Independent,  B  and  W  being  one  dimensional  Wiener  processes 

and  Xq  Is  real.  An  approximate  computation  procedure  (at  periodic  sampling 

points  t  *  nh,  n  c  IN ,  h  >  0)  of  the  unnormalized  conditional  measure  ot  of 

the  Zakai  equation  Is  proposed.  It  Is  shown  that  E  l<Wf>  -  *  0  « 

-1/2 

n  *  for  t  *  nh  and  the  speed  of  convergence  Is  0(n  '  ).  (Here  f  Is  a 
sufficiently  regular  function  and  o'  Is  the  solution  of  the  approximating 
recursive  equation). 
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HUI-HSIUNG  KUO 


1.  Applications  of  differential  geometry  In  statistical  problems  of  stochastic 
processes  [1]  (jointly  with  G.  Kalllanpur). 

This  work  Is  motivated  by  the  recent  paper  of  Amari  [Ann,  of  Statist.  8, 

1980]  on  the  use  of  differential  geometric  techniques  In  (finite-dimensional) 
parameter  estimation  In  statistics. 

With  the  Idea  of  extending  the  elegent  approach  of  Amari  and  Efron  to 
stochastic  processes,  we  study  the  concepts  of  Riemmanian  connection,  Christoffel 
symbols,  curvature  and  geodesics,  "curved"  families  of  measures,  etc.,  for 
Rlemnanlan  Hilbert  manifolds.  The  reason  for  considering  Hilbert  manifolds 
Is  that  the  natural  parameters  In  Inference  In  stochastic  processes  are  Infinite 
dimensional  manifolds.  A  preliminary  investigation  has  been  made,  largely 
consisting  of  gathering  known  facts  on  Hilbert  manifolds.  We  have  written  up 
a  unified  account  of  that  part  of  Amari 's  work  that  Is  dimension  free  and  have 
tried  to  formulate  the  statistical  problems  In  stochastic  processes  to  which 
these  concepts  are  to  be  applied.  The  part  already  done  Is  being  prepared  as 
a  Technical  Report.  Further  work  will  be  carried  out  later. 

2.  On  generalized  Brownian  functionals  [2]  (jointly  with  G.  Kalllanpur) 

The  regularity  property  of  Donsker's  delta  function  6(B(t)-x)  Is  established. 
The  result  Is  applied  to  obtain  a  simple  proof  of  Ito's  lemma  for  f (B(t) )  where 
f  is  a  tempered  distribution. 
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HANNU  0.  NIEMI 


2 

Prediction  properties  of  stationary  random  fields  on  Z  [1]. 

o 

Let  X  ,  (m,n)  e  Z  be  a  second  order  stationary  random  field  with 
in  t  n 

spectral  measure  The  work  Is  devoted  to  establishing  the  relationship 
between  the  spectral  representation  of  the  random  field  and  the  four-fold 
Wold  decomposition  for  Xm  n  obtained  by  Kallianpur  and  Mandrekar  (Prediction 
Theory  and  Harmonic  Analysis,  V.  Mandrekar  and  H.  Salehl,  eds..  North  Holland, 
1983,  165-190.)  The  differences  between  the  two  Wold-Halmos  decompositions 
(Theorems  1.1  and  1.2)  of  the  cited  paper  have  been  expressed  in  spectral  terms 
A  spectral  characterization  of  the  horizontal  (resp.  vertical)  Wold  decompo¬ 
sition  for  XJn  n  has  also  been  obtained,  clarifying  earlier  results  by 
Chiang-Tse  Pel  (Theor.  Probability  Appl .  2,  1957). 
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MOHSEN  POURAHMADI 

Or.  Pourahmadl  continued  his  work  on  prediction  of  multivariate  stationary 
processes  and  also  worked  on  prediction  problems  for  stable  processes  and  on 
nonlinear  time  series  models. 

1.  On  Minimality  and  Interpolation  of  Harmonlzable  Stable  Processes  [1] 

It  Is  shown  that  a  harmonlzable  symmetric  a-stable  process,  l.e.  the 

Fourier  coefficients  of  a  process  with  Independent  symmetric  a-stable  Increments, 

with  spectral  density  w  Is  minimal.  If  and  only  if  w>0  a.e.  and  w-1/(a-l)  Is 
Integrable.  Algorithms  for  the  linear  Interpolator  and  Interpolation  error  of 

such  processes  are  given  when  several  values  of  the  process  are  missing.  The 

algorithm  for  the  linear  Interpolator  Is  convergent  If  w  satisfies  Huchenhoupt's 

(Ao)  condition. 

2.  On  Statlonarlty  of  the  Solution  of  a  Doubly  Stochastic  Model  [2]. 

Consider  the  discrete  parameter  process  {X^>  satisfying  the  doubly 
stochastic  model  Xt  ■  +  et,  where  {<Pt>  and  {c^}  are  also  stochastic 

processes.  Necessary  and  sufficient  conditions  on  are  given  for  {Xt>  to  be 
a  second  order  process.  When  Is  a  strictly  stationary  process,  some  sufficient 
conditions  in  terms  of  are  given  which  guarantee  the  wide  sense  stationarity 
of  <Xt).  It  turns  out  that  for  these  problems  the  distribution  and  dependence 
structure  of  the  process  {log|4>t|}  play  an  Important  role. 

3.  On  the  Rate  of  Wean  Convergence  of  Finite  Linear  Predictors  of  Multivariate 
Stationary  Stochastic  Processes  [3] 

Consider  a  multivariate  weakly  stationary  stochastic  process  {Xn)  with  the 
spectral  density  matrix  f  satisfying  the  boundedness  condition.  It  is  shown  that 
If  the  entries  of  f  are  analytic  functions  of  0  on  [-u,tt],  then  the  rates  of 


convergence  of  the  one-step  ahead  linear  least  squares  predictor  of  (Xn>  based 
on  a  finite  segment  of  the  past,  and  of  the  partial  sum  of  the  Infintle 
linear  least  squares  predictor  of  the  process  to  the  Kolmogorov-Wiener  predictor, 
are  at  least  exponential. 

4.  Infinite  Order  Autoregressive  Representations  of  Multivariate  Stationary 
Stochastic  Processes  [4] 

Consider  a  q-varlate  weakly  stationary  stochastic  process  (Xn)  with  the 
spectral  density  W.  The  problem  of  autoregressive  representation  of  (Xn)  or 
equivalently  the  autogreresslve  representation  of  the  linear  least  squares 
predictor  of  XR  based  on  the  Infinite  past  Is  studied.  It  is  shown  that  for 
every  W  In  a  large  class  of  densities,  the  corresponding  process  has  a  mean 
convergent  autoregressive  representation.  This  class  Includes  as  special  sub¬ 
classes,  the  densities  studied  by  Masanl  (1960)  and  the  author  (1984).  As  a 
consequence  It  Is  shown  that  the  condition  of  Integrablllty  of  W'1  or  minimality 
of  (Xn)  Is  dispensable  for  this  problem.  When  W  Is  not  In  this  class  or  when  W 
has  zeros  of  order  2  or  more.  It  Is  shown  that  In  this  case  (Xn>  has  a  mean  Abel 
suimtable  or  mean  compounded  Cesaro  summable  autoregressive  representation. 
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HOLGER  ROOTZEN 


Dr.  Rootzen  completed  his  research  with  J.  Sternby  on  a  Bayesian  approach 
to  consistency  In  least  squares  estimation  (cf.  abstract  In  1983  report).  This 
work  has  now  been  published  ([1]). 

In  the  area  of  extreme  values.  Dr.  Rootzen  considered  the  greatest  possible 
rates  of  convergence  of  the  (normalized)  maximum  to  an  extreme  value  distribu¬ 
tion,  for  1.1. d.  sequences.  He  showed  that  the  difference  between  this  distribu¬ 
tion  of  the  maximum  and  Its  limit,  could  not  tend  to  zero  at  "arbitrary 
exponential  rates",  unless  the  distribution  of  the  terms  of  the  sequence  Is 
Itself  of  extreme  value  type.  However,  examples  show  that  any  specific  exponential 
rate  Is  attainable  by  some  sequence  having  a  non-extreme  value  distribution. 

This  work  has  been  reported  In  the  publication  [2]. 
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A  symmetric  stable  process.  X(t),  t  e  T,  can  be  represented  as  a  stochastic 
integral  /J  f(t,s)dM(s),  t  «  T,  where  M  is  a  stable  motion,  and  a  fundamental 
problem  is  to  relate  properties  of  the  process  X  with  properties  of  the  kernel 
f.  Here,  we  assume  that  the  sample  paths  of  the  process  X  belong  to  a  certain 
Banach  space  V(T)  of  functions  on  T.  Vie  show  that  the  kernel  f  can  be  modified 
in  such  a  way  that  Its  sections  f(‘,s)  belong  to  the  same  Banach  space  V(T) 

(the  converse  Is  false  in  most  interesting  cases).  This  enables  us  to  use  some 
techniques  of  probabilities  on  Banach  spaces  to  estimate  the  moments  of  the 
norm  of  the  sample  paths  of  the  process: | | X( • ) J 1 y(T) *  Applying  this  result  we 
solve  the  problem  of  estimation  of  the  moments  of  a  double  integral  of  stable 
motion.  He  also  characterize  Infinitely  divisible  distributions  whose  tail 
behaviour  Is  similar  to  the  tail  behaviour  of  stable  distributions.  This 
permits  us  to  extend  previous  results  to  a  more  general  class  of  distributions, 
including  in  particular  all  semistable  distributions;  and  also  to  establish  the 
equivalence  of  the  moments  of  processes  with  stochastic  integral  representations 
In  terms  of  motions  determined  by  these  more  general  distributions,  with  the 
corresponding  moments  of  stable  processes. 


2.  Moment  inequalities  for  p-stable  stochastic  Integrals  [2]. 


f  p 

popents  of  the  random  scale  change  T(t,w)  *  |F(s,oi)|  ds.  We  also  generalize 

and  apply  a  result  of  Bass  and  Cranston  (Ann.  Probability,  11,  1983,  578-588) 
concerning  moments  of  exit  times  and  exit  points  for  symmetric  stable  processes. 

3.  Continuity  of  some  random  Integral  mappings  [3]. 

Jurek  and  Vervaat  (Z.  Wahr.,  62,  1983,  247-262)  proved  that  every 

self decomposable  measure  y  on  a  Banach  space  B  Is  the  distribution  of  the 

stochastic  integral  e“tdX(t),  where  the  Banach  space  valued  process  X  has 

stationary  Independent  Increments  and  E  log+l|X(l)||  <  ~.  Moreover, 

the  mapping  which  takes  the  law  of  X(1)  to  the  law  of  /q  e~ldX(t)  is  an  algebraic 

Isomorphism  of  the  convolution  semigroup  ID^0g(B)  of  all  Infinitely  divisible  y's 

with  /log+| | x | |dy(x)  <  ®  and  the  Levy  class  L(B).  In  this  paper  we  solve  the 

problem  of  the  continuity  of  this  Isomorphism.  We  find  a  metric  topology  of 

IDlog(B)  such  that  this  mapping  becomes  a  topological  Isomorphism.  This  enables  us 

to  determine  generators  of  the  Levy  class  L(B).  These  results  are  also  generalized 

to  the  class  of  operator  self decomposable  measures.  We  also  study,  as  a  necessary 

tool,  the  relationship  between  the  uniform  integrability  of  Infinitely  divisible 

* 

probability  measures  and  of  their  Levy  measures,  which  may  be  of  independent  interest. 
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DAG  TJ0STHEIM 

Dr.  Tjdstheim  continued  his  research  in  the  area  of  nonlinear  models 


for  time  series. 

1.  Some  Doubly  Stochastic  Time  Series  Models  [1] 

We  consider  time  series  models  obtained  by  replacing  the  parameters  of 
autoregressive  models  by  stochastic  processes.  Special  attention  is  given  to 
the  problem  of  finding  conditions  for  stationarity  and  to  the  problem  of 
forecasting.  For  the  first  problem  we  are  only  able  to  obtain  solutions  in 
special  cases,  and  the  emphasis  is  on  techniques  rather  than  obtaining  the  most 
general  results  In  each  case.  For  the  second  problem  more  complete  results  are 
obtained  by  exploiting  similarities  with  discrete  time  (nonlinear)  filtering 
theory.  The  methods  introduced  are  illustrated  on  two  standard  examples,  one 
of  state  space  type  and  one  where  the  parameter  process  is  a  Markov  chain. 

2.  Recent  Developments  in  Nonlinear  Time  Series  Modelling  [2]. 

We  review  developments  In  nonlinear  time  series  that  have  taken  place  in 
the  last  decade  or  so.  Three  main  model  classes  are  considered:  Models 
motivated  by  nonlinear  differential  equations,  bilinear  models  and  doubly 
stochastic  models.  For  these  classes  we  study  problems  of  existence  of  a 
stationary  solution,  Invertlbillty,  prediction,  model  fitting  and  quite  briefly 
statistical  estimation.  We  also  discuss  some  attempts  of  constructing  a  more 
general  framework  for  nonlinear  time  series  analysis. 

3.  least  Squares  Estimates  and  Order  Determination  Procedures  for  Autoregressive 
Processes  with  a  Time  Dependent  Variance  (with  J.  Paulsen)  [3]. 

We  study  nonstationary  autoregressive  processes,  where  the  variance  of 
the  generating  white  noise  process  Is  allowed  to  depend  on  time.  It  Is  shown 
that  ordinary  least  squares  estimates  are  strongly  consistent  and  with  a 


proper  scaling  factor  asymptotically  normal,  but,  as  can  be  expected,  they  are 
not  efficient.  Furthermore,  AIC  type  order  determination  criteria,  used  as  If 
the  underlying  process  Is  stationary,  are  consistent,  whereas  Identification 
of  order  In  terms  of  the  partial  autocorrelation  function  may  lead  one  astray. 

4.  Estimation  In  Nonlinear  Time  Series  Models  I:  Stationary  Series  [4], 

A  general  framework  for  analyzing  estimates  In  nonlinear  time  series  models 
Is  developed.  Ergodlc  strictly  stationary  series  are  treated.  General  conditions 
for  strong  consistency  and  asymptotic  normality  are  derived  both  for  condi¬ 
tional  least  squares  and  maximum  likelihood  type  estimates.  Examples  are 
taken  from  exponential  autoregressive,  random  coefficient  autoregressive  and 
bilinear  time  series  models.  Some  nonstationary  models  and  examples  are 
treated  in  a  sequel  to  this  paper. 

5.  Estimation  In  Nonlinear  Time  Series  Models  II:  Some  Nonstationary  Series  [5] 

In  [4]  a  general  framework  was  Introduced  for  analyzing  estimates  in 
stationary  nonlinear  time  series  models.  Here  the  framework  Is  enlarged  to 
Include  certain  nonstationary  and  nonlinear  series.  General  conditions  for 
strong  consistency  and  asymptotic  normality  are  derived  both  for  conditional 
least  squares  and  maximum  likelihood  type  estimates.  Examples  are  taken  from 
threshold  autoregressive,  random  coefficient  autoregressive  and  doubly  stochastic 
(dynamic  state  space)  models.  The  emphasis  in  the  examples  Is  on  conditional 
least  squares  estimates. 
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RESEARCH  IN  STATISTICAL  ESTIMATION  AND  INFERENCE 


RAYMOND  J.  CARROLL 


During  the  past  year,  research  has  continued  In  the  areas  of  transformations, 
heteroscedastlclty  and  related  statistical  topics.  Leonard  A.  Stefanski  and 
David  M.  Glltlnan  completed  their  Ph.D.  dissertations  In  December,  1983,  under 
the  direction  of  Professor  Carroll.  Mr.  Douglas  Simpson  and  Ms.  Marie  Davldian 
are  working  on  their  Ph.D.  dissertations  under  the  direction  of  Professor 
Carroll;  Simpson  should  finish  In  June  1985,  while  Davldian  has  just  started 
her  work  and  should  finish  In  June  1986. 

1.  Influence  and  Measurement  Error  In  Logistic  Regression.  The  Ph.D.  thesis  of 
Leonard  A.  Stefanski,  supervised  by  Professor  Carroll. 

This  dissertation  concerns  the  use  of  logistic  regression  when  certain 
standard  model  assumptions  are  violated.  Chapters  I  and  II  study  the  problem 
of  estimating  regression  parameters  when  covariates  are  subject  to  measurement 
error.  The  latter  chapters  study  robust  methods  applicable  to  logistic 
regression. 

To  facilitate  study  of  the  errors-ln-variables  problem  a  small  measure¬ 
ment  error  asymptotic  theory  is  developed.  This  allows  comparison  of  certain 
estimators  which  have  appeared  In  the  literature  and  also  suggests  new  estimators 
which  are  shown  to  have  better  asymptotic  properties.  A  small  Monte-Carlo 
study  confirms  the  superiority  of  the  new  estimators  In  certain  settings.  In 
the  course  of  studying  the  asymptotic  behavior  of  the  various  estimators 
Interesting  use  Is  made  of  some  random  convex  analysis. 

To  deal  with  the  problem  of  messy  data,  l.e.  outliers  and  extreme  co¬ 
variables,  several  bounded  Influence  estimators  are  proposed.  The  optimality 
properties  of  these  estimators  are  studied  In  Chapter  III.  Asymptotic  theory 
for  the  robust  procedures  Is  given  In  Chapter  IV.  Finally,  Chapter  V  concludes 
the  thesis  with  an  application  of  these  methods  to  two  sets  of  data. 


43 


2.  Transformations  In  Regression:  A  Robust  Analysis  (with  David  Ruppert). 

We  consider  two  approaches  to  robust  estimation  for  the  Box-Cox  power 
transformation  model.  One  approach  maximizes  weighted,  modified  likeli¬ 
hoods.  A  second  approach  bounds  the  self-standardized  gross-error  sensitivity, 
a  measure  of  the  potential  Influence  of  outliers  pioneered  by  Krasker  and 
Welsch  (JASA,  1982). 

Among  our  primary  concerns  Is  the  performance  of  these  estimators  on 
actual  data.  In  examples  that  we  study,  there  seem  to  be  only  minor  differences 
between  these  three  estimators,  but  they  behave  rather  differently  than  the 
maximum  likelihood  estimator  or  estimators  that  bound  only  the  influence  of 
the  residuals. 

Confidence  limits  for  the  transformation  parameter  can  be  obtained  by 
using  a  large-sample  normal  approximation  or  by  modified  likelihood-ratio 
testing. 

These  examples  show  that  model  selection,  determination  of  the  transforma¬ 
tion  parameter,  and  outlier  Identification  are  fundamentally  Interconnected. 

3.  Bounding  Influence  and  Leverage  In  Logistic  Regression  (with  L.A.  Stefanski 
and  D.  Ruppert). 

Logistic  regression  Is  very  sensitive  to  discordant  data.  Resistant 
fitting  procedures,  though  clearly  desirable,  have  not  received  their  due 
attention.  This  paper  studies  two  bounded  Influence  estimators  similar  In 
spirit  to  the  Krasker-Wel sch  estimator  for  linear  regression.  In  motivating 
the  bounded  Influence  procedures  some  theoretical  optimality  results  for 
Krasker-Wel sch  type  estimators  Are  obtained.  These  results  are  not  specific 
to  logistic  regression  and  generalize  to  other  statistical  models.  In  particular 
the  class  of  generalized  linear  models.  The  proposed  estimation  procedures 
are  Illustrated  with  examples. 
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4.  Bounded  Influence  Estimation  In  Heteroscedastic  linear  Models.  The  Ph.D. 
thesis  of  David  M.  Giltlnan,  supervised  by  Professor  Carroll. 

In  a  heteroscedastic  linear  model.  If  there  Is  no  replication,  the  usual 
approach  Is  to  model  the  variance.  A  common  situation  Is  that  where  the 
variance  Is  modelled  as  a  function  of  the  mean  response,  or  of  a  subset  of  the 
explanatory  variables.  Maximum  likelihood  estimation  of  the  variance  parameter 
In  such  models  Is  very  sensitive  to  extreme  data  points.  In  this  dissertation 
alternative  'bounded  Influence1  methods  are  developed  which  overcome  this 
problem. 

The  model  for  the  variance  which  Is  considered  Is  given  by  the  relationship 
o.j  *  exp[h'(T.j)eJ,  where  h  Is  a  known  vector-valued  function  and  8  Is  the  unknown 
variance  parameter.  Using  this  model,  there  Is  a  parallel  between  the  problem 
of  developing  efficient  bounded  Influence  estimators  for  £  and  that  of  developing 
efficient  bounded  Influence  estimators  for  the  regression  parameter  In  the 
homoscedastlc  regression  case.  Extension  of  several  such  methods  of  estimating 
6  to  Include  the  problem  of  estimating  6_  Is  discussed. 

The  method  proposed  by  Krasker  (1980)  Is  generalized  to  Include  estima¬ 
tion  of  8^  and  shown  to  be  optimal  In  a  certain  sense.  The  question  of  existence 
of  a  solution  to  the  estimating  equations  Is  discussed. 

The  methods  of  Krasker  and  Welsch  (1982)  In  the  homoscedastlc  regression 
case  are  extended  to  cover  estimation  of  e.  A  three-stage  weighted  regression 
estimate  based  on  these  techniques  Is  proposed  and  the  associated  Influence 
calculations  are  presented. 

The  class  of  bounded  Influence  regression  estimators  proposed  by  Mallows 
(1975)  Is  considered.  A  necessary  condition  for  a  strongly  optimal  weight 
function  is  derived  and  the  resulting  estimator  of  £  is  generalized  to  obtain 
an  estimator  of  8.  Mallows  estimators  are  shown  to  possess  a  stability  of 


variance  which  Is  absent  In  Krasker-Welsch  estimators.  A  three-stage  Mallows 
regression  estimate  Is  proposed  and  Its  Influence  function  Is  given. 


Computation  of  the  three-stage  estimates  Is  discussed,  and  the  performance 
of  these  methods  Is  evaluated  by  using  them  In  the  analysis  of  a  number  of 
data  sets. 

5.  A  Note  on  Levene's  Tests  for  Equality  of  Variances  (with  H.  Schneider). 

Consider  testing  for  equality  of  variance  In  a  one-way  analysis  of  variance. 
Levene's  test  Is  the  usual  F-test  for  equality  of  mean  computed  on  pseudo¬ 
observations,  which  one  defines  as  the  absolute  deviations  of  the  data  points 
from  an  estimate  of  the  group  "center".  We  show  that,  asymptotically,  Levene's 
test  has  the  correct  level  whenever  the  estimate  of  group  "center"  is  an 
estimate  of  group  median.  This  explains  why  published  Monte-Carlo  studies 
have  found  that  Levene's  original  proposal  of  centering  at  the  sample  mean  has 
the  correct  level  only  for  symmetric  distributions,  while  centering  at  the 
sample  median  has  correct  level  even  for  asymmetric  distributions.  General¬ 
izations  are  discussed. 

6.  Some  New  Estimation  Methods  For  Weighted  Regression  When  There  Are  Possible 
Outliers  (with  D.M.  Giltlnan  and  D.  Ruppert) 

We  consider  the  problem  of  fitting  a  linear  regression  when  the  data  have 
nonconstant  variances,  l.e.  the  problem  of  heteroscedastic  regression.  We  choose 
as  a  model  that  the  variances  depend  on  a  function  of  known  predictor  variables 
and  unknown  parameters.  We  propose  two  new  estimators  which  have  bounded 
Influence.  These  estimators  are  based  on  the  Idea  of  (1)  separately  bounding 
all  sources  of  Influence,  the  Mallows  approach,and  (2)  bounding  all  sources  of 
Influence  simultaneously,  the  Hampel  and  Krasker/Welsch  approach.  We  find 
that  bounding  the  effect  on  predicted  values  gives  better  estimators  in  practice 


46 


than  other  methods.  The  new  estimators  are  evaluated  on  sets  of  data  which 
contain  outliers. 


7.  Some  Diagnostics  for  Accelerated  Life  Testing  (with  H.  Schneider  and 
L.  Nelssfeld) 

In  accelerated  life  testing,  we  are  often  concerned  with  understanding  the 
failure  distribution  of  a  component  at  a  fixed  level  of  a  predictor  variable  such 
as  temperature,  this  level  being  the  one  which  will  be  encountered  in  practice. 

In  many  Instances  It  Is  Impractical  to  do  experimentation  at  the  level  or 
temperature  of  Interest,  because  of  time  constraints.  What  Is  often  done  in 
this  case  Is  to  relate,  through  a  model,  failure  of  a  component  to  the  predictor 
variable.  We  have  encountered  situations  within  the  context  of  accelerated 
life  testing  In  which  one  or  more  components  seem  to  fall  much  too  early.  In 
this  paper  we  describe  some  simple  diagnostics  for  determining  whether  single 
or  pairs  of  points  have  an  enormous  impact  on  the  Inferences  at  the  level  of 
the  predictor  variable  which  Is  of  real  interest. 
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HELMUT  SCHNEIDER 

Dr.  Schneider  visited  during  the  1983-84  academic  year  and  conducted 
research  In  the  following  areas: 

1.  A  numerical  analysis  of  sampling  plans  based  on  prior  distribution  and  costs 
This  paper  deals  with  some  pitfalls  linked  with  the  sampling  model  based 
on  prior  distribution  and  costs.  First,  a  model  Is  designed  which  encompasses 
most  of  the  existing  Bayesian  cost  models.  The  efficiency  of  sampling  plans 
Is  Investigated  In  a  numerical  study.  It  Is  shown  that  under  realistic  assump¬ 
tions,  described  by  Dodge  (1969)  and  Schilling  (1982),  sampling  plans  based 
on  prior  distributions  and  costs  are  only  efficient  In  an  outlier  model, 

1. e.  If  almost  all  lots  are  of  good  quality  and  only  a  low  number  of  lots, 
denoted  as  outlier  lots,  have  very  poor  quality.  Furthermore,  It  Is  demon¬ 
strated  that  for  the  Polya  distribution  a  gain  of  sampling  Is  linked  with  a 
high  percentage  of  rejections,  l.e.  when  the  prior  distribution  cost  relation¬ 
ship  Is  such  that  less  than  5*  of  the  lots  should  be  rejected, sampling  be¬ 
comes  Inefficient. 

2.  The  Performance  of  Variable  Sampling  Plans  when  the  Normal  Distribution  Is 
Truncated. 

The  robustness  of  standard  variable  sampling  plans  by  Lleberman  and 
Resnlkoff  Is  considered  with  respect  to  a  truncation  of  the  normal  distribution. 
It  Is  shown  how  variable  sampling  plans  can  be  designed  If  the  truncation  point 
and  o2  are  known. 

3.  Analysis  of  a  Multi-Item  Inventory  Problem  Using  Optimal  Policy  Surfaces 

An  Interactive  algorithm  Is  presented  which  allows  the  selection  of  a 
combination  of  service-level.  Investment,  workload  and  storage  room  which 
Is  "appropriate"  for  the  management.  This  method  produces  combinations  which 
lie  on  an  optimal  surface. 


Three  methods  for  applying  the  EM  algorithm  to  censored  data  are  considered, 
the  Buckley- James  (1979),  a  proposed  simpler  nonparametric  method  and  a 
normal  model  for  censored  data.  A  new  estimator  for  the  variance  of  y  in  the 
Buckley- James  model  Is  proposed  and  simulations  comparing  the  three  methods  are 
described.  To  Illustrate  the  use  of  these  methods  they  are  applied  to  the 
Stanford  heart  transplant  data. 
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diagnostics  at  the  annual  meeting  of  the  American  Statistical  Association. 
Throughout  the  year  he  served  as  an  associate  editor  for  the  Journal  of  the 
American  Statistical  Association  and  the  Annals  of  Statistics,  as  well  as  a 
member  of  a  National  Academy  of  Sciences  panel  on  Youth  Employment  in  the 
United  States. 

G.  Kalllanpur  gave  an  Invited  lecture  at  the  I.F.I.P.  Conference  on  Stochastic 
Differential  Systems  &  Filtering  at  Lumlny-Marseilles  France,  (March  84),  and 
a  series  of  talks  on  nonlinear  filtering  and  stochastic  differential  equations  at 
the  Department  of  Mathematics,  University  of  Alberta,  Edmonton,  Alberta,  Canada 
(April  84).  He  also  served  as  editor  of  Sankhya  and  as  associate  editor  of  the 
Journal  of  Multivariate  Analysis,  the  Journal  of  Applied  Mathematics  and 
Optimization,  and  Stochastic  Processes  and  Their  Applications. 

M.R.  Leadbetter  gave  Invited  talks  at  the  Oberwolfach  conference  on  order 
statistics  and  quantiles  (April  84)  and  at  the  annual  Conference  on  Stochastic 
Processes  and  Their  Applications  held  at  Gothenberg,  Sweden  (June  84).  He 
continued  his  activity  on  the  Council  for  the  Institute  of  Mathematical 
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Statistics,  and  represented  the  IMS  on  the  American  Mathematical  Society 
commltte  for  surnner  conferences. 

T.  Funakl  gave  colloquium  talks  at  Rutgers  University  and  New  York 
University. 

W.  Hudson  gave  a  talk  at  the  Conference  on  Probability  on  Banach  Spaces 
held  at  Tufts  University  In  Boston  (July  84). 

R.  Karandlkar  gave  colloquium  talks  at  Michigan  State  University  and  at 
the  Operations  Research  Department  of  the  University  of  North  Carolina. 

M.  Pourahmadl  attended  the  NSF  Conference  on  Inference  for  Stochastic 
Processes  held  at  Lexington,  KY  (May  84). 

H.  Rootzen  gave  colloqula  at  Colorado  State  University  and  at  Johns 
Hopkins  University. 

J.  Ros Inski  gave  a  talk  at  the  Conference  on  Probability  on  Banach  Spaces 
held  at  Tufts  University  in  Boston  (July  84). 

H.  Schneider  presented  seminars  at  McMaster  University,  the  University 
of  Calgary  and  the  University  of  North  Carolina  at  Chapel  Hill. 

D.  Tjdsthelm  gave  a  talk  at  the  Princeton  Conference  on  Information 
Sciences  and  Systems  held  at  Princeton  University  (March  84). 

A.  Heron  gave  a  talk  at  the  Conference  on  Probability  on  Banach  Spaces 
held  at  Tufts  University  In  Boston  (July  84). 


